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Abstract 



We study the nature of the two-dimensional quantum critical point separating 
two phases with and without long-range spin-density- wave order, which has 
been recently observed in cuprate superconductors. We consider the Landau- 
Ginzburg- Wilson Hamiltonian associated with the spin-density critical modes, 
perform a mean-field analysis of the phase diagram, and study the correspond- 
ing renormalization-group flow in two different perturbative schemes at five 
and six loops, respectively. The analysis supports the existence of a stable 
fixed point in the full theory whose basin of attraction includes systems with 
collinear spin-density-wave order, as observed in experiments. The stable fixed 
point is characterized by an enlarged 0(4) <g> 0(3) symmetry. The continuous 
transition observed in experiments is expected to belong to this universality 
class. The corresponding critical exponents are v = 0.9(2) and rj = 0.15(10). 

PACS: 05.10.Cc, 75.30.Kz, 74.72.-h, 05.70.Jk. 
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I. INTRODUCTION 



In the last few decades several aspects of cuprate superconductors (SCs) have been 
studied and many efforts have been spent to understand the unique and complex phase 
diagram exhibited by this class of materials; see, e.g., Ref. 1. Superconductivity in cuprates 
appears to be due to a mechanism analogous to the BCS one in ordinary superconductors. 
However, superconductivity is only one of the characteristic features of these materials. 
There are many other new properties that require more complex mechanisms and can be 
understood only if the interplay between BCS and additional order parameters is considered. 
For instance, at T ps 0, L^-^Sr^CuC^ at very low doping 5 is an insulator with long-range 
magnetic order. Increasing 5, at 5 ~ 0.055 an insulator-superconductor first-order transition 
takes place, giving rise to a superconducting state in which spins are still magnetically 
ordered 2 . At 5 ~ 0.14 another phase transition occurs, and, for 5 > 0.14, the material shows 
no magnetic order — it is paramagnetic — but is still superconducting. Neutron-scattering 
experiments 3 suggested that this transition is continuous. Moreover, in the ordered phase 
5 < 0.14, they revealed the presence of collinearly polarized spin-density waves (SDWs) with 
wavevectors 



where 9 is a function of the doping concentration and a is the lattice spacing. The wave 
vectors Kj are two-dimensional since cuprates are supposed to be made of weakly interacting 
planes and thus behave approximately as two-dimensional systems. Following Ref. 4, we 
assume that superconductivity is not relevant at the transition which is instead driven by 
the interaction among the SDW degrees of freedom. Since T ps one should take into 
account the quantum nature of the system. Quantum phase transitions can be studied by 
introducing a supplementary dimension parametrized by an imaginary time variable r. The 
relevant order parameter is the spin field which is parametrized as 



where $ a j are complex amplitudes. There are two interesting limiting cases. The first 
one is when the order parameter can be written as $ a (r, r) = e lola n a , which corresponds 
to collinearly polarized SDWs. The second one is when $ a (r, r) = n 0j i +i n a2 , with 
n a ,i ■ n aj2 = and |n 0j i| = |n 0j2 |, which corresponds to circularly polarized SDWs. In 
cuprates experiments indicate that the ground state shows a collinear behavior. 5 

The standard strategy for writing down an effective Hamiltonian for a given physical 
system consists in considering all polynomials of the order parameter of order less than 
or equal to four that are compatible with the expected symmetries. In the SDW-SC-to- 
SC phase transition the order parameter is the complex field 3> j(r, t), with a = 1,2 and 
% = 1,2,3. The corresponding symmetries are the following: (i) SO(3) spin rotations: 
^ai Oij^ a j; (ii) Translational symmetry of the spin waves: $ a j — > e %OLa $ a j; (iii) Spatial 
inversion: $ a j — > (iv) Interchange of the 1 and 2 axes: $ij <-> $2« and x <-> y. The most 
general Hamiltonian with these symmetries is 4 
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n = J d 2 rdT{\d T <b l \ 2 + v 2 l \d x <$> l \ 2 + vl\d y <b 1 \ 2 + |<M> 2 | 2 + 

+v 2 \d x <£> 2 \ 2 + vf\d y & 2 \ 2 + r(|$!| 2 + |$ 2 | 2 ) + 
+^(l$i| 4 +l$ 2 | 4 ) + ^(|^| 2 + |^| 2 ) + 

+wi i0 |$i| 2 |$ 2 | 2 + w 2i0 |$i • $ 2 | 2 + w 3fi \®l ■ $2| 2 } , (1.3) 

where v 1 and v 2 are parameters called SDW velocities. Terms such as $* -d T ^ a are forbidden 
by spatial inversion symmetry and terms like i$*-d x $ a , even if permitted by all symmetries, 
can be eliminated by redefining the fields as $ a — ► e* qa ' r $ a . Hamiltonian (1.3) admits several 
different ground states depending on the values of the parameters. They are classified in 
App. A. In particular, there is the possibility that both fields correspond to collinearly 
polarized SDWs as observed in experiments: $1 = e iai ni and $ 2 = e ia2 n 2 , where the 
vectors ri! and n 2 satisfy either n x = n 2 or • n 2 = 0. 

In this paper we investigate the nature of the fixed points (FPs) of the renormalization- 
group (RG) flow of the effective Hamiltonian (1.3). If a stable FP exists and its attraction 
domain includes systems with collinearly polarized SDWs, then the SDW-SC-to-SC tran- 
sition may be continuous. Otherwise, it must be of first order. In our study, we consider 
only the case v 1 = v 2 that simplifies the analysis and allows us to perform a high-order 
perturbative analysis. Therefore, we consider the theory 

H= j d d x (Id^l 2 + \d^ 2 \ 2 ) + r(|^| 2 + |$ 2 | 2 ) + ^(l^l 4 + |$ 2 | 4 ) 

+ ^(l$ 2 | 2 + l$2| 2 ) + «>l,0|*l| 2 |*2| 2 + «*,0|*1 ' $ 2| 2 + W 3fi m • $ 2 | 2 | , (1.4) 

where the field <3> a j is a complex 2 x N matrix, a — 1, 2, % — 1, . . . , N. The physically relevant 
case is = 3. 

We first perform a standard analysis close to four dimensions, 6 computing the RG func- 
tions in powers of e = 4 — d. A one-loop analysis indicates that a stable FP exists only 
for N > 42.8. Apparently, this result casts doubts on the existence of a stable FP in three 
dimensions. However, in three dimensions there may exist FPs that are absent for e< 1. 
This is indeed what happens in the Ginzburg-Landau model of superconductors, in which 
a complex scalar field couples to a gauge field 7 and in 0(2)®0(n) symmetric models. 8 ' 9 
Thus, a more careful investigation of the RG flow in three dimensions calls for strictly three- 
dimensional perturbative schemes. For this purpose we consider two field-theoretical pertur- 
bative approaches: the minimal-subtraction scheme without e expansion 10 (in the following 
we will indicate it as 3c?-MS scheme) and the massive zero-momentum (MZM) renormal- 
ization scheme. 11 The use of two different schemes is crucial, since the comparison of the 
corresponding results provides a nontrivial check on the reliability of our conclusions. In the 
3<i-MS scheme one considers the massless (critical) theory in dimensional regularization, 12 
determines the RG functions from the divergences appearing in the perturbative expansion 
of the correlation functions, and finally sets e = 4 — d = 1 without expanding in powers of 
e (this scheme therefore differs from the standard e expansion 6 ). In the MZM scheme one 
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considers instead the three-dimensional massive theory in the disordered (high-temperature) 
phase. We compute the f3 functions to five loops in the 3rf-MS scheme and to six loops in 
the MZM scheme. We use a symbolic manipulation program that generates the diagrams 
(approximately one thousand at six loops) and computes their symmetry and group factors, 
and the compilation of Feynman integrals of Refs. 13, 14. The series are available on request. 
The perturbative expansions are then resummed using the known large-order behavior. 

The perturbative analysis of the RG flow in the full theory is not sufficiently stable to 
provide reliable results. Therefore, we have focused on the stability of the FPs that occur 
in specific submodels of Hamiltonian (1.4). The analysis of the perturbative series indicates 
the stability of the 0(4) <g> 0(3) collinear FP that occurs in the model with w 1:0 = u lj0 — -u 2 ,o 
and u>2,o = ^3,0 = ^2,0 < 0. Moreover, its basin of attraction includes systems with collinear 
SDWs. Therefore, we expect the continuous transition observed experimentally in cuprates 
to belong to this universality class. This implies an effective enlargement of the symmetry 
at the transition point. The corresponding critical exponents would be 

v = 0.9(2), = 0.15(10). (1.5) 

The paper is organized as follows. In Sec. II we discuss the possible ordered phases that 
occur in model (1.3) in the mean-field approximation. Details are given in App. A. In Sec. Ill 
we discuss the FP structure close to four dimensions in the standard e expansion. Sec. IV 
contains the main results of this work. We consider three different submodels (Sec. IV A) 
and then investigate the stability properties of the FPs occuring in each of them (Sections 
IV B, IV C, and IV D). Conclusions are presented in Sec. V. In App. B and C we give some 
technical details. 

II. MEAN-FIELD ANALYSIS 

The phase diagram of Hamiltonian (1.3) can be studied in the mean-field approximation. 
Due to the large number of couplings the analysis is quite complex. We have limited our 
considerations to the case N < 3. We summarize here the results that are derived in App. A. 
For r > the system is disordered and $1 = $ 2 = 0. For r = a continuous phase transition 
occurs followed by a magnetized phase with r < 0. The nature of the ordered phase depends 
on the values of the quartic parameters. The analysis reported in App. A shows that there 
are seven possibilities: 

(1) $1 is a collinear SDW ($1 = e ioil n, n real) while $ 2 = 0. 

(2) $1 is a circularly polarized SDW (<3>i = e mi (ni + in 2 ), ni and n 2 real, |ni| = |n 2 |, 
rix ■ n 2 = 0) while $ 2 = 0. 

(3) $1 and $ 2 correspond to collinear SDWs with the same axis and amplitude: $1 = 
e iai n, $ 2 = e ia2 n, n real. 

(4) $1 and $ 2 correspond to collinear SDWs with orthogonal axes and same amplitude: 
$1 = e tai n.i, $ 2 = e* a2 n 2 , ni,n 2 real, ni ■ n 2 = 0, |ni| = |n 2 |. 
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(5) $i and $2 are circularly polarized SDWs with the same rotation plane and amplitude: 
$ 1 = e iai (ni + in 2 ), $ 2 = e M2 (rii + in 2 ), ni and n 2 real, |ni] = |n 2 |, ■ n 2 = 0. 

(6) $1 is a collinear SDW and $ 2 is a circularly polarized SDW. The rotation plane of $ 2 
is orthogonal to the axis of $1. Explicitly: $1 = e tai ni, $ 2 = e ia2 (n 2 + m 3 ), rij real, 
|n 2 | = |n 3 |, rij • rij = 0. 

(7) $1 and $ 2 are elliptically polarized SDWs with different rotation planes but with the 
same amplitude, |$i| = |$ 2 |. 

For cuprates the relevant solutions are (3) and (4). Necessary conditions to obtain (3) are 
w 2 ,o + w 3fi < and 

wi,o + w 2 ,o + w 3fi - u lfi < u 2fi < Min [u lfi - w lj0 - w 2fi - w 3fi , -w 2 ,o, -w 3fi ], (2.1) 

while (4) requires w 2t o + w 3 $ > and 

wi,o ~ ui,o < U2,o < Min [u lj0 - w lj0 , w 2fi , w 3fi ]- (2.2) 

These conditions are not sufficient, since for some values of the parameters satisfying 
Eqs. (2.1) or (2.2) the ordered phase is given by solutions (6) or (7). Note that the sign 
of u 2>0 is not the relevant parameter that selects the collinear SDWs among all possible 
solutions. 

It is interesting to note that the mean-field solution predicts either $i|| ( I > 2 or $i±$ 2 in 
the case of collinear SDWs. This result is easy to understand. If both fields correspond 
to collinear SDWs, then one can take $1 and $ 2 real. In this case the only term of the 
Hamiltonian that contains a scalar product of the two fields is (u> 2 ,o + w 3t o)(&i ■ $ 2 ) 2 that 
forces the two fields to be either parallel or orthogonal, depending on the sign of w 2t o + w 3j o- 
Note that this also holds if we add additional higher-order terms to the Hamiltonian, as long 
as the transition is continuous. Indeed, for a continuous transition $ a — > at the transition 
($ a = in the disordered phase) and thus higher-order terms do not play any role. On the 
other hand, this relation may not be valid if the transition is of first order. Also the coupling 
to the charge- density waves (CDWs) that are present in cuprates 15 ' 4 does not change this 
conclusion, since they couple to the scalars |$ a | 2 - 

Solutions (3) and (4) also satisfy |$i| = |$ 2 |. This property does not necessarily hold if 
we take into account the CDWs (see Refs. 15, 4 for an extensive discussion). Indeed, let <p\ 
and 2 be the complex amplitudes of the CDWs coupled respectively to $ 2 and & 2 . In the 
absence of the CDW-SDW coupling, for some values of the CDW Hamiltonian parameters, 
the ordered solution corresponds to |0i| 7^ 0, <p 2 = 0. If now the CDW-SDW coupling is 
included, one may obtain a ground state with |0i| 7^ \<p 2 \ 7^ and |$i| 7^ |$ 2 | 7^ 0. 

III. RG FLOW CLOSE TO FOUR DIMENSIONS 

The RG flow close to four dimensions can be investigated perturbatively in e = 4 — d. 
In the minimal-subtraction (MS) the one-loop (3 functions are: 
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P U1 = -eui + (N + A)u\ + Au x u 2 + <±u\ + Nw\ + w\ + w\ + 2wiw 2 + 2w x w 3 , 
At 2 = ~ eM 2 + QuiU 2 + Nu\ + 2w 2 w 3) 

f3 Wl = —ewi + 2w\ + w\ + W3 + 2(iV + l)-uiu?i + 4w 2 wi + 2u\w 2 + 2u±w 3: 
(3 W2 = —ew 2 + Nw% + 2u 1 w 2 + Au 2 w 3 + 4w 1 w 2 + 2w 2 w 3) 

(3 W3 = -ew 3 + Nwl + 2uiw 3 + Au 2 w 2 + Awiw 3 + 2w 2 w 3 , (3.1) 

where Ui, Wi are the renormalized quartic couplings corresponding to the quartic Hamiltonian 
parameters w^o, ^i,o- They are normalized so that, at tree level, g = gofi~ e / Ad, where g and g 
label the renormalized and Hamiltonian parameters respectively and = 2 d ~ 1 it d l 2 T(d/2). 
The FPs of the RG flow are the common zeroes of the (3 functions. For N = 3 there are 
4 FPs while for N = 2 there are 7 FPs: they are all unstable. Only for N > 42.8 does a 
stable FP exist. It has u 2 = w 2 = w 3 (for iV — > 00 we obtain u\ = u 2 = w 2 = w 3 = e/N, 
Wi — 0), so that at the FP the symmetry becomes 0(4) <g> O(N). This FP is the chiral FP 
that occurs in O(M) ® O(N) in the large- N limit. 16 

In order to determine the behavior in three dimensions, one should extend the compu- 
tation to higher order in e and determine the function N c (e) = 42.8 + 0(e) such that the 
chiral FP point identified above exists for iV > iV c (e) and is no longer present for smaller 
values of N. We have not pursued this approach for several reasons. First, the analogous 
five-loop computation that was performed in the O(N) ® 0(2) model 17 ' 16 ' 18 was not able to 
explain the correct physics of these models for iV = 2, 3 (see Sec. II. D in Ref. 8). Moreover, 
this calculation is only concerned with the stable FP that is present for e = (in the present 
case the chiral 0(4) ®0(iV) FP), while in d = 3 the stable FP may be different, an unstable 
or even a new FP. The analysis that will be presented in the next Section favors this last 
possibility. 



IV. SUBMODELS AND THEIR STABILITY 

The three-dimensional properties of the RG flow are determined by its FPs. Some of them 
can be identified by considering particular cases in which some of the quartic parameters 
vanish. The corresponding FPs are also FPs of the general theory. In this section, we identify 
some of them, and then determine their stability with respect to the complete theory. 



A. Some particular cases 

For particular values of the couplings Hamiltonian (1.4) reduces to that of simpler models. 
Three cases have already been extensively studied in the literature: 19 

(1) For wi t o = w 2 fl = w 3t Q = there is no interaction between the two SDWs and Hamilto- 
nian (1.4) reduces to that of two identical decoupled 0(2)(g>0(iV)-symmetric models. 
The general 0(m) <g> 0(n)-symmetric model is defined by the Hamiltonian density 17,19 



+ rc 
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02,0 

4! 



'm'raj 



where a j is a real n x m matrix field (a = 1, 
(4.1) is obtained from Eq. (1.4) by setting $ ai 

«1,0 = ^1,0/3 - 02,o/6, «2,0 = 02,0/6, Wi l0 = W 2 , = ^3,0 = . 



, n and % = 1, . . . , m) 
^ + and 



(4.1) 



Hamiltonian 



(4.2) 



The properties of 0(2) ® O(iV) models are reviewed in Refs. 8, 17, 19, 20. In three 
dimensions perturbative calculations within the MZM scheme 21,22 and within the 3d- 
MS scheme 8 indicate the presence of a stable chiral FP with attraction domain in the 
region g 2 ,o > for all values of N (only for N = 6 the evidence is less clear, since the 
MZM analysis does not apparently support it). For N = 2, these conclusions have 
been recently confirmed by a Monte Carlo simulation. 8 A stable collinear FP for g 2 < 
exists for N < 4. 9 ' 23 Apart from the collinear FP for N = 2, these FPs do not exist 
close to four dimensions. For N = 2 the collinear FP is equivalent to an XY FP and 
corresponds to g* = g* X y> 9% = ~9*xyi where g* XY is the FP value of the renormalized 
coupling in the 0(2) 4 model. 

(2) For wifi = Uifi — w 2 ,o an d w 2 ,o = ^3,0 — u 2,o, Hamiltonian (1.4) reduces to (4.1) with 
m = 4 and n = N. The correspondence is given by 



li 



$2i = 



where 



V2 ' 
is a 4 x N matrix, and 

£1,0 = 3(111,0 + 112,0), #2,0 



V2 



Qu 



2,0 • 



(4.3) 



(4.4) 



We have already discussed the FPs of the 0(4) ® 0(2) theory. The 0(4) ® 0(3) theory 
does not present stable FPs for gi > 0. 24 Analyses of the available six-loop series in 
the MZM scheme and five-loop series in the 3-d MS scheme indicate the presence of a 
stable collinear FP for g 2 < 0. 25 This FP does not exist close to four dimensions. 



(3) For U2fl = u>2,o = ^3,0 — we obtain the van model with n - 
so-called mn model is defined by the Hamiltonian density 26 ' 19 



2 and m = 2N. The 



5^(<9 M 0ai) 2 + r(j 



+ 



92,0 

4! 



5> 



at] 



2 .6 2 - 



(4.5) 



where <fi a i is a real nxm matrix, i.e., a 
is obtained by setting 

$ai = {<t>ai + i(t>a,i+N)/V2, 



1, . . . , n and i — 1, . . . , m. The correspondence 



01,0 = 3lUi n, 02,0 = 3(«i )0 - Win) 



(4.6) 



A stable FP is the 0(m) FP with g\ = and g 2 = <?o(m)' wnere 9o{m) ls ^he FP value of 
the renormalized coupling in the 0(m)-symmetric vector model. In three dimensions 
the analysis of five- and six-loop series 27 indicates the presence of a second stable FP 
with g 2 < for n = 2 and m — 2, 3, and 4. 
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Beside these three models, there are two other submodels for which no results are avail- 
able: 

(a) For w 2 fi = w 3;0 — we obtain two chiral models coupled by an energy-energy term. 
Note that in this model the RG flow does not cross the planes u 2 = and w 1 = 0. 

(b) For w 2 fi = ^3,o — w o we obtain a model with an additional U(l) symmetry: $ x — > <3> x , 
$2 — ► $2- m this m odel the RG flow does not cross the plane w — 0. 

Finally, note an additional symmetry of Hamiltonian (1.4). It is invariant under $ x — > $*, 
$ 2 - ► ^2) w 2,o — * w 3,o an d w 3,o — * w 2,o, while the other couplings are unchanged. This 
implies that the RG flow in the space of renormalized couplings does not cross the plane 
w 2 = w 3 and that, for any FP with w 2 > w 3 there is an equivalent one with w 2 < w 3 . In 
particular, we can limit our considerations to w 2 > w 3 . 

In order to study the RG flow of the theory one can start by discussing the stability in 
the full theory of the FPs of the models (1), (2), and (3) discussed above. 

For N = 2 and N = 3, the only cases we consider, model (1) has two FPs: 

(la) the chiral FP, in which g 1 = #* ch and g 2 = g* 2 ch ; correspondingly u\ = 0j )Ch /3 - ^ )Ch /6, 
u *2 = #2,ch/ 6 > °> w l = w^ = w* 3 = 0; 

(lb) the collinear FP, in which gi = g^ cl and g 2 = cl ; correspondingly u\ = g\ cl /3— g 2 c i/6, 
u *2 = 9tJ Q < °> w l = w 2 : = w 3 * = 0. 

Here gi^ and gi iC \ are the chiral and collinear FPs of the 0(2)<g)0(iV) theory. The analogous 
FPs are present in model (2): 

(2a) the chiral FP, in which g 1 = gl ch and g 2 = g 2ch ; correspondingly u\ = <? lch /3 — g 2ch /6, 
u* 2 = g 2>ch /6 < 0, wl = ul - ul, wl = W3 = u* 2 ; It does not exist for N = 3. This is the 
FP that is relevant for iV > N c (e) rs 42.8 + 0(e) close to four dimensions; 

(2b) the collinear FP, in which g\ = (?* cl and g 2 = g^^; correspondingly u\ = g\ cl /3— g\ cl /6, 
u 2 = g 2jd /6 < 0, wl = ul - u 2 , w 2 = w 3 = u 2 . It exists for both N = 2 and N = 3. 

Here g^ and g^i are the chiral and collinear FPs of the 0(4)®0(iV) theory. Finally, the 
mn theory gives two FPs: 

(3a) the 0(2N) FP. This is unstable in the full theory, being already unstable in model (1); 

(3b) the mn FP g x = g ljmn , g 2 = g 2:fnn ; correspondingly u\ = g* tmn /3 + 3 2 * mn /3, wl = 
9i,mn/3, u 2 = w 2 = = 0. It exists only for N = 2. 

In the following we study the stability of these FPs in the complete theory (1.4). For 
this purpose, using the (5 functions of the general theory we have computed the stability 
matrices of the FPs at six and five loops respectively in the MZM and 3<i-MS schemes. The 
perturbative series have been resummed by using the conformal-mapping method described, 
e.g., in Ref. 28. For a FP belonging to a submodel, the large-order behavior needed for the 
conformal-mapping summation is the same as that characterizing all series of that submodel. 
For all submodels we consider, the large-order behavior is already known. 21 > 24 > 29 
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B. Stability of the decoupled 0(2) <g> 0(N) fixed points 



We want to establish the stability properties of the decoupled 0(2) <g> 0(N) FPs (la) 
and (lb) in the complete theory (1.4). For this purpose we need the RG dimensions of the 
operators present in Hamiltonian (1.4) that break the symmetry of model (1), i.e., of the 
operators associated with the quartic couplings Wi. It is useful to rewrite them as 

wi, |$i| 2 |$ 2 | 2 + u>2,o|$i ' $2| 2 + w 3fi \<5>l ■ $ 2 | 2 = W 00 P 00 + W U P U + WoA, (4.7) 
where 

Woo = w lfi + ^(^2,o + w 3fl ), W n = -w 2 ,o + w 3j0 , W 02 = w 2 ,o + w 3j0 , (4.8) 

p _ n (lW2) p n (2) p n (2) ,, Q N 

^00 = <^00 ^00 > Ml = U ll,ij U ll,iji ^02 = 2_^ U 02,ij U 02,ij> l 4 - y J 

ij 

and, using the correspondence (4.2), 

Off = E * { M\ 0% = - (4.10) 

ei 

ifl) - y 6 {a) 6 {a) - — 5 V <£ {a V a) 



N 

ek 



ek 



The quadratic operator O^J transform as a spin-m and a spin-/ representation with respect 
to the 0(2) and O(N) groups, respectively. Since P o, Pu, and P 02 belong to different 
irreducible representations, they do not mix under RG transformations at the decoupled 
0(2) ® O(N) FPs. Their RG dimensions Y nd can be derived from the RG dimensions y m i of 
the quadratic operators at the 0(2) <g> O(N) FP, using the relation 

Y ml = 2y ml -3. (4.11) 

The quadratic term Off corresponds to the energy operator and thus yoo = 1/^ and Y 00 = 
a/u, where a and v are the specific-heat and correlation-length critical exponents of the 
given 0(2) <g> O(N) FP. The RG dimensions y n and y 02 were computed in Ref. 23 (there, 
they are named y\ and y 3 respectively). 
At the chiral FP (la) we obtain: 30 

Y 00 = 0.3(3) Y u = 1.6(3) Y 02 = 0.04(8) for N = 3; 

Yoo = 0.2(3) Yn = 1.9(4) Y 02 = -0.4(2) for N = 2. 

At the collinear FP (lb) we obtain: 30 

Yoo = 0.3(2) Y n = -0.6(2) Y 02 = 1.0(3) for N = 3; 

Yoo = -0.2182(8) Y n = -0.022(8) Y 02 = 0.9240(11) for N = 2. 

These results show that the decoupled 0(2) ®0(N) FPs are unstable in the complete theory 
(1.4) for both N = 3,2. 

It is also interesting to discuss submodels (a) and (b) mentioned in Sec. IV A. In model 
(a) one should only consider P o- The numerical results apparently indicate that the FPs 
are always unstable (but, with the present errors, we cannot really exclude the opposite 
possibility), except in one case. For N = 2, the collinear FP is stable. In model (b) one 
should consider P o and P 2- For N = 2, 3, all FPs are unstable. 
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FIG. 1. Distribution of the results for Y\ (left) and Yz (right) obtained by varying the resum- 
mation parameters a and b as a function of the number of loops in the MZM and 3d-MS schemes. 
Here iV = 3. 



C. Stability of the 0(4) ® 0(N) FPs 

Here we investigate the stability of FPs (2a) (it does not exist for N = 3) and (2b). 
For this purpose we must compute the RG dimensions of the perturbations of the 0(4) (g) 
O(N) model appearing in the complete theory. This is done in App. B. There are two 
relevant operators with RG dimensions Y\ and Y 2 . The corresponding perturbative series 
are reported in App. B. They are analyzed using the conformal mapping method. 28,31 The 
errors we will report takes into account the variation of the estimates when changing the 
resummation parameters b, a denned in Ref. 31 — we use b = 3, . . . , 18 and a = 0, . . . , 4 — and 
the uncertainty of the FP coordinates. 

The analysis of the sixdoop series in the MZM scheme and of the five-loop 3d-MS series 
gives the following results at the collinear FP (2b): 25 
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Y 1 = 


-0.4(4) 


Y 2 = 


-0.95(7) 


for 


N = 


2 


(MZM), 


Y 1 = 


-0.6(9) 


Y 2 = 


-1.2(1.0) 


for 


N = 


2 


(3d-MS), 


Y 1 = 


-1.5(1.2) 


Y 2 = 


-0.42(10) 


for 


N = 


3 


(MZM), 


Yx = 


-0.8(1.5) 


Y 2 = 


-0.1(2) 


for 


N = 


3 


(3d-MS). 



(4.12) 



For N = 3 the MZM and 3rf-MS results are consistent and apparently indicate that Y\ 
and Y 2 are negative, though with somewhat large errors. A better understanding of the 
relevance of the two operators can be obtained from Fig. 1, where we give the distributions 
of the estimates of Y 1 and Y 2 obtained by varying the parameters a and b. For Yi low-order 
calculations predict Y ± > 0. However, as the number of loops increases, Y x decreases. The 
six-loop MZM results indicate that Y x < 0, a result that is also supported by the trend 
observed in the 3<i-MS results. As for Y 2 , the MZM results clearly indicate Y 2 < 0. However, 
this is not fully confirmed by the 3d-MS results. Though they give Y 2 < 0, there is a trend 
towards larger values of Y 2 . Overall, these results support the stability of the collinear FP 
(2b) in the complete theory for N = 3. 

Similar conclusions hold for N = 2. For completeness, we report here the corresponding 
critical exponents: 33 



V = 


0.71(7) 


7] 


= 0.12(1) 


for 


N = 


2 


(MZM); 


V = 


0.76(10) 


V 


= 0.11(6) 


for 


N = 


2 


(3d-MS); 


V = 


0.89(16) 


V 


= 0.18(3) 


for 


N = 


3 


(MZM); 


V = 


0.88(22) 


V 


= 0.10(10) 


for 


N = 


3 


(3d-MS) 



(4.13) 



For iV = 2 we also study the stability of the chiral FP (2a). Using the results of Refs. 8, 22 
for the FP we have: 32 

Yl = -0.03(7) Y 2 = 0.9(2) (MZM); 

Y 1 = -0.2(3) Y 2 = 0.73(15) (3d-MS). 1 ' ' 

The chiral FP is clearly unstable. 

Finally, note that the same discussion also applies to submodel (b), since the stability of 
the FP depends on the same operators with RG dimensions Y 1 and Y 2 . For submodel (a) 
one should only consider Y\. In this case also the chiral FP (2b) might be stable. 



D. Stability of the mn FP for N = 2 

Here we investigate the stability of FP (3b) for N — 2 (it does not exist for N — 3). 
For this purpose we must compute the RG dimensions of the perturbations of the mn FP 
appearing in the complete theory. This is done in App. C. There are two relevant operators 
with RG dimensions Yi and Y 2 . The analysis of the perturbative series in the MZM scheme 
gives 

Y 1 = -4.0(2.6), Y 2 = -0.6(2). (4.15) 

The results in the 3<i-MS scheme are very imprecise, although negative values for Y\ and Y 2 
seem to be favored. There results indicate, although with limited confidence, that the mn 
FP present for iV = 2 may be stable in the complete theory. 
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V. CONCLUSIONS 



In this paper we have studied the quantum phase transition that occurs in two- 
dimensional systems that exhibit an ordered phase with SDW order. The effective Hamilto- 
nian of the relevant critical modes & ai is given in Eq. (1.3). A detailed mean- field analysis 
shows that in some parameter region Hamiltonian (1.3) has a continuous transition separat- 
ing a spin disordered phase from an ordered phase characterized by two collinearly polarized 
SDWs. There are two different possibilities for the the axes of these SDWs: either <&i||$2 or 
$i_L$2- We have then investigated the role of fluctuations in a simplified model in which the 
two SDWs have the same velocity. For this purpose we have generated six-loop perturbative 
series in the MZM scheme and five-loop series in dimensional regular izat ion with minimal 
subtraction. Close to four dimensions, an analytic e-expansion calculation shows no presence 
of stable FPs. However, past experience indicates that a FP may exist in three dimensions 
and be absent for e <C I. Therefore, we have considered two strictly three-dimensional 
schemes. We have analyzed the stability of some FPs that belong to known submodels. The 
analysis of the perturbative series supports the stability of the 0(4)® 0(3) collinear FP. The 
analyses of the MZM and 3d-MS expansions do not provide sufficiently stable results for the 
RG flow in the full theory, i.e. in the general space of its five quartic couplings. In particular, 
they do not allow us to draw any definite conclusion on the existence of other stable FPs. 
In any case, even without the analysis of the full flow, simple considerations (reported in 
App. A) show that systems with collinear SDWs with the same axis (the mean-field solution 
(3) reported in Sec. II) are in the attraction domain of the 0(4) £g> 0(3) collinear FP. 

It should be remarked that our RG analysis is only valid for v\ = v 2 - In order to extend 
the results to the generic case v\ ^ v 2 one should also consider the operator 

O v = \d X ^\ 2 - Idy^l 2 - \d X <S> 2 \ 2 + |^<I> 2 | 2 (5.1) 

and determine its RG dimension y v at the 0(4) <g> 0(3) collinear FP. If y v < the previous 
conclusions are unchanged. On the other hand, if y v > the 0(4) <g> 0(3) collinear FP is 
unstable with respect to the perturbation O v . In this case the transition may be of first 
order or continuous depending on the existence and attraction domain of a stable FP with 
Vi 7^ v 2 . Note that, from a practical point of view, our results are of interest even if y v > 0. 
Indeed, one expects the SDW velocities v 1 and v 2 to be close in magnitude, of the order of 
the spin-wave velocity of the Neel state of the undoped insulator. 4 Therefore, the RG flow 
always starts very close to the stable FP of the theory with v\ — v 2 , and thus the critical 
behavior is controlled by this FP except in a narrow interval around the critical doping. 

Experiments indicate that the SDW-SC-to-SC transition is continuous and is associated 
with collinear SDWs. It is thus natural to conjecture that its critical behavior is controlled 
by the 0(4) <g> 0(3) collinear FP, since this FP is stable in model (1.4) and its basin of 
attraction includes systems with collinear SDWs. The corresponding critical exponents are 
then predicted to be v = 0.9(2), rj = 0.15(10). 
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APPENDIX A: GROUND-STATE CONFIGURATIONS 

In this Appendix we compute the possible ground-state configurations of Hamiltonian 
(1.4), that allow us to identify the possible symmetry-breaking patterns. We consider 
translation-invariant configurations and the space-independent Hamiltonian density 

if ($1, $ 2 ) = r(|$!| 2 + |$ 2 | 2 ) + if 4 ($i, $2), (Al) 

where H4 is the part of the Hamiltonian that is quartic in the fields. Since H4 > for 
stability, for r > the ground state always corresponds to $1 = $ 2 = 0. For r < 0, 
$1 = $ 2 = is a local maximum of H and thus the ground state is nontrivial. The value 
r = corresponds to a second-order transition point in the mean-field approximation. In 
order to determine the ground states for r < 0, we will first determine all stationary points 
of H; the ground state is the one with the lowest energy. Note that, if $1, $ 2 is a stationary 
point, then 

$ 2 ) = ^(|$ 1 | 2 + |$ 2 | 2 ) = -/J 4 ($ 1 ,$ 2 ). (A2) 
This relation is quite general. Indeed, assume H to be of the form 

H = ^ njfafa + ^2 9ijki<fii<fij<fik<fii- (A3) 

ij ijkl 

Then 

H = - + 2 r n^i = 2 fagi: ~ 9ijki<t>i4>j4>k(t>i- (A4) 

i ij i ijkl 

On a stationary solution, the derivative vanishes, proving Eq. (A2). 

The calculation of the ground states also allows us to determine the stability domain 
of the Hamiltonian. Indeed, a point in the coupling space does not belong to the stability 
domain if there is a field such that H4 < 0. Being H4 homogeneous, it is not restrictive to 
consider only fields such that |$i| 2 + |$ 2 | 2 = 1- Thus, the determination of the minima of 
H 4 is equivalent to the determination of the minima of (H — r) where r is now interpreted 
as a Lagrange multiplier. Eq. (A2) shows that H 4 can be negative only for r > 0. Thus, the 
stability domain of H4 is obtained by determining the stationary points of H for r positive. 

In order to determine the minima, we can use the symmetry of the Hamiltonian. Using 
the 0(N) symmetry we can always write 

Re$i = (a,0, . . .), Im$i = (6, c, 0, ...)■ (A5) 
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Then, by using the U(l) symmetry we can also fix b = 0. Indeed, we first perform an 0(2) 
rotation on the first two components: 

ip[ = ipi cos 9 + ip 2 sin 9 ip' 2 = — ipi sin# + ip 2 cos 9, (A6) 

where ip is either Re^i or Im$!. Then, we apply a U(l) rotation, $'/ = e"*^. If we choose 

26c a sin 9 , . „. 

tan2#= — t ana = — -, A7 

cr + ¥ — c z o sin — c cos 9 

the transformed field has the form (A5) with 6 = 0. Once $i has been fixed we can use 
0(N — 2) and U(l) rotations to write 

$ 2 = (d + ie, g + if,l + ih, im, 0, . . .)• (A8) 

If A < 3, one can use U(l) rotations to set I = 0. 

The analysis of the minima is nontrivial due to the complexity of the stationarity equa- 
tions. We have only consider the case A < 3 that is relevant experimentally. Other ground 
states are present for A > 4. We found seven relevant minima (only five of them occur for 
A = 2): 

1) a 2 = -r/u 12 , H = -r 2 /(2u 12 ). 

2) a 2 = c 2 = -r/ui, H = -r 2 /(2 Ul ). 

3) a 2 = d 2 = -r/(u 12 + w 1 + w+), H = -r 2 /(u l2 + w 1 + w+). 

4) a 2 = j 2 = -r/(«ia + w 1 ),H= -r 2 /{u 12 + Wl ). 

5a,b) a 2 = c 2 = d 2 = f 2 = -r/(2u 1 + 2wi + w + ± w-), H = -2r 2 j(2u x + 2w x + w+ ± w-). 

6) [A > 3] c 2 = -r(ui - Wl )/A e , d 2 = h 2 = -r{u l2 - w 1 )/{2A 6 ), H = -r 2 ( Ul + u 12 - 
2w 1 )/(2A 6 ), A 6 = mu 12 - w\. 

7) [A >3,w±^ 0] a 2 = d 2 ^ 0, c,f,h^ 0, ad/(cf) = w + /w_, a 2 + c 2 = d 2 + f 2 + h 2 , with 
energy H = -r 2 [u 2 (w 2 + w 3 ) + w 2 w 3 }/ A 7 , with A 7 = (u 12 + w 1 )w 2 w 3 + w + u 2 (ui +w 1 ). 
Alternatively, if we define the four vectors ti = Re^i, t 2 = Im<E>i, t 3 = Re $2; 
^4 = Im$ 2 , and tij = ti ■ tj, the solution can be characterized more geometrically as 
follows: tn = t 22 = t 33 = tu = H/ (2r), t 12 = t 3A , £13 = t 2A , and t 14 = t 23 , with 

j.2 222/A2 j.2 222/A2 j.2 2 2 2/a2 /in\ 

t 12 = r w 2 w 3 /A 7 t 13 = r u 2 w 2 /A 7 t 14 = r u 2 w 3 /A 7 . (A9) 

Whenever a component is not explicitly written, it vanishes. Moreover, we defined u\ 2 = 
u\ + -u 2 , w + = w 2 + w 3 , W- = w 2 — w 3 and we simplified the notation writing u\ instead of 
1*1,0, etc. Beside the seven solutions reported above, for w- ^ we also found stationary 
points with e = g = h = and 

df -4u| -w 2 _+w 2 + ± a/( 4m 2 + w 2 _ - w\) 2 - \<ou\w 2 _ 
ac Au 2 W- 
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A numerical analysis indicates that they are never absolute minima of the Hamiltonian 
and thus they are never relevant for the ground-state calculation. For this reason, these 
solutions have not been included above. The computation of all stationary points is quite 
straightforward, except for solution 7. We shall now briefly sketch how it is derived. Assume 
that e = g = and a, c, d, /, h ^ and define E a = (l/a)dH/da, etc. Then, 

f 

J -(E f - E h ) = -adw- + cfw + = 0. (All) 

Using this relation, we can rewrite E a , E c , E d , and Ef as linear equations in a 2 , c 2 , d 2 , f 2 , 
and h 2 . Considering also a 2 d 2 w 2 _ = c 2 f 2 w 2 + that follows from Eq. (All), we obtain a system 
of equations that allows us to determine all components. 

Given the list of solutions, we can determine the stability domain of the Hamiltonian. 
Using solutions 1-5, we obtain the necessary conditions 

Ui > 0, U12 > 0, -Ui2 + U>l>0, 

U12 + w\ + w + > 0, ui + w\ + 2^- w + ^ u '-) > 0- (A12) 

These conditions are sufficient for iV = 2. For iV > 3 we must also consider solutions 6 and 
7. Solution 6 gives the necessary condition 

wi > -y/uiUi 2 . (A13) 

Numerically, we find that solution 7 is also relevant for stability, although we have not been 
able to write down an easy condition. 

For cuprates the relevant solutions are 3 and 4. In view of the possibility that the 0(4) (8) 
0(3) FP is stable it is important to understand to which ground state of the 0(4) ® 0(3) 
Hamiltonian (4.1) they correspond. For generic n and m, m > n, model (4.1) is stable for 
gifi > and ngip — (n — 1)^2,0 > and has two ground states depending on the sign of g2,o : 
f° r 92,o > the ground state is chiral, while for g 2 ,o < the ground state is collinear. 17 The 
corresponding energies are: 

H = r (chiral); 

21*101,0 - (n - 1)52,0] 

(A14) 

3r 2 

H = (collinear). 

2<7i,o 

Using U12 = 51,0/3, u 2 = 52,0/6, wi = (51,0 - 52,o)/3, w + = g 2 ,o/3, W- = 0, we immediately 
see that for m = 4 and n = 2 and n = 3 solutions 1 and 3 correspond to the collinear 
case. For n = 3 solutions 6 and 7 correspond to the chiral case, while solutions 2, 4, and 
5 correspond to a stationary state that is never a ground state in the chiral theory. For 
n = 2 instead, solutions 2, 4, and 5 are those corresponding to the chiral case. This result is 
relevant to identify the attraction domain of the 0(4) <g> 0(3) collinear FP in the full theory. 
Indeed, it shows that the attraction domain of this FP includes systems whose ground state 
is given by solutions 1 and 3 (and therefore two collinearly polarized SDWs). Nothing can 
be said on the other solutions: in this case an analysis of the RG flow of the full theory is 
needed. 
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APPENDIX B: RENORMALIZATION-GROUP DIMENSIONS OF THE 
PERTURBATIONS AT THE 0(4)<g>0(7V) FIXED POINTS 



We need to classify the operators that break 

0(4) © 0(N) -> (U(l) © U(l)) © 0{N) (S0(2) © S0(2)) © 0(iV) 



(Bl) 



This is essentially discussed in Ref. 23. There are, however, two differences: first, we have 
only SO (2) symmetry, instead of 0(2) symmetry; second, there is an additional exchange 
symmetry that forbids the appearance of spin-2 operators. In the notations of Ref. 23 (M 
and iV of Ref. 23 correspond to N and 4 respectively) we define 



Pi - oftS + o 

P<1 



(4,4) 



1144 + ^2233 + ^2244 J 
/-0(4,r) . n (4,r) ^(4,r) n (4,r) 
L^iqiq "I" ^1414 ~T~ ^2323 ^2424 5 



,(4,4) 



(B2) 



'1313 



-^3 = ^1234 5 



where e j is the real field defined in Eq. (4.3). Note that P 3 would be forbidden if we had 
0(2) invariance instead of SO (2) invariance. Moreover, P 2 and P 3 correspond to different 
components of the same operator, so that they have the same RG dimension. 
Hamiltonian (1.4) can then be written as 



n = J ddx E \ K v ^) 2 + + *i(E ^ 2 + ** E 

a a a,b 

+t 3 P 1 +t 4 P 2 + t 5 P 3 , 



(B3) 



where 



fi 

*5 



— (2-Ui + 2-u 2 + wi + u> 2 + u> 3 ) 
1 

36 
1 

12 
1 



(ui + Au 2 — wi + 2w 2 + 2w 3 ) 
'—Ui - u 2 + W 1 + w 2 + w 3 ) 



— (-2u 1 + 4u 2 + 2w l 
^{w 3 - w 2 ) . 



w 2 - w 3 ) 



(B4) 



Since all operators are irreducible with respect to 0(4) © 0(N) transformations, if the 
couplings belong to the 0(4) ©O(iV) theory, the stability matrix defined with respect to the 
couplings U has the form 



/fin Q 21 \ 

Q21 fi 2 2 

fii 

fi 2 

\ tt 2 / 



(B5) 



16 



Here Y\ = — fii is the RG dimension of O^tfi an d ^2 = —^2 is the RG dimension of 0^ r \ 
In the MZM scheme for iV = 2 we find: 

Q ± = 1 - (0.238732 Ml + 0.31831 u 2 ) + 

+ (0.0324838 uj + 0.0614494 Ul u 2 + 0.0342428 m^) + 

— (0.00570145 m? + 0.016556 m?m 2 + 0.0218393 u x u\ + 0.00451785 m^) + 

+ (0.00156398 u\ + 0.0052904 m?m 2 + 0.0097641 u\u\ + 0.00663506 u x u\ + 

+0.000911563 m*) - (0.00045218m? + 0.00188296 m?m 2 + 0.00448727 m? m^ 

+0.00487778 m? u 2 + 0.00209129 u x u\ + 0.000250855 m^) + (0.000158666 m? + 

+0.00074771 lu\u 2 + 0.00212573 u\u\ + 0.00311938 m? Ml) + 0.00230887 m? M^j + 

+0.000759722 m^ + 0.0000778366 m^), (B6) 

Q 2 = 1 - (0.238732 ui - 0.159155 u 2 ) + 

+ (0.0324838 m? - 0.0145415 u x u 2 - 0.0164178 m^) + 

-(0.00570145 m? - 0.00426028 u\u 2 - 0.00222856 u x u\ + 0.00294388 u\) + 
+ (0.00156398 u\ - 0.000683968 u\ u 2 + 9.94721 • 10" 6 uj u\ + 0.000950042 Ul ul + 
— 0.000152361 m 2 ) - (0.00045218m? - 0.000214037m? u 2 + 0.00011481 u\u\ + 
+0.000428198 u\ u 2 - 0.0000869867m 1 m 2 1 - 0.0000326009 ul) + (0.000158666 m? + 
-0.000038529 m?m 2 + 0.000128979 m?m^ + 0.000233286 m?m 2 ; - 0.0000225715 m?m^ + 
-0.0000164528 u x u\ + 8.9658 • 10~ 6 u 6 2 ). (B7) 

For N = 3 we obtain: 

fti = 1 - (0.238732 Mi + 0.397887m 2 ) + 

+ (0.0378782 m? + 0.0687202 Ul u 2 + 0.0413963 u\) + 

-(0.00623158 m? + 0.0193945 uj u 2 + 0.0290399 u lU 2 2 + 0.00434652 ul) + 

+ (0.0020362m? + 0.00594275 m?m 2 + 0.0126485 uju 2 2 + 0.00940724 Ul ul + 

+0.000356974 u%) - (0.000572885 m? + 0.00220348 m?m 2 + 0.00583754 m?m^ 

+0.0069323 m?m 2 : + 0.00263869 u\ u 2 + 0.000106569 m^) + (0.000227195 m? + 

+0.000881969 m?m 2 + 0.00277192 m?m^ + 0.00439833 m? ul + 0.00335526 uju 4 2 + 

+0.000841582 Ul u\ + 0.0000352387 m^), (B8) 

tt 2 = l- (0.238732 U\ - 0.0795775 m 2 ) + 

+ (0.0378782 u\ - 0.00727073 u x u 2 - 0.0314283 u\) + 

-(0.00623158 m? - 0.00255553 m?m 2 - 0.00424767 u x u\ + 0.00445854 ul) + 

+ (0.0020362 m? - 0.000329962 m?m 2 - 0.000353201 m?m 2 : + 0.00115746 u x ul + 

-0.000132244 u\) - (0.000572885m? - 0.000151507m? m 2 + 0.000153088 m?m^ + 

+0.000587068 m? m 2 - 0.0000150878 m u\ - 0.0000901012 ul) + (0.000227195 m? + 

-0.0000200211 m?m 2 + 0.000205186 m?m^ + 0.000321211 m?m^ + 0.000031431 uj u\ + 

-0.0000663864 Mi ul + 0.0000193321 u 6 2 ). (B9) 
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In the 3d-MS we find for N = 2: 

Q 1 = 1- (6m 1 + 8m 2 ) + 

+ (30.5 u\ + 58 mi m 2 + 32 u\) + 

-(327.297m? + 928.74 m?m 2 + 1185.79^^ + 258.397 m^) + 
+ (5835.31 u\ + 20132.4 u\ u 2 + 35648.3 u\ u\ + 23296.2 m u\ + 
+4377.86 u\) - (123668 m? + 506531 u\u 2 + 1.1389- 10 6 m?m^ + 
+1.1852- 10 6 ^^ + 552355mi m^ + 85949 u|), (BIO) 

Q 2 = l- (6«i-4u 2 ) + 

+ (30.5 m? _ 14 Mi m 2 - 16 u|) + 

-(327.297m? - 219.798m? u 2 - 163.596mim^ + 103.301 u\) + 

+ (5835.31 u\ - 3001.7m? m 2 - 2298.63 m^ + 1415.85 m^ + 

+ 198.484 u\) - (123668 m? - 48887.2 m?m 2 - 24466.2 m?m^ + 

+34690.3 u\ u\ - 12396.4 Ul u\ - 4256.31 u\). (Bl 1) 

and finally for N = 3: 

ft x = 1 - (6mi + 10m 2 ) + 

+ (35.5 m? + 65 Mi m 2 + 38.5 uj) + 

— (369.646 m? + 1082.64 m? m 2 + 1548.49 ii\ u\ + 251.598 m 2 ) + 

+ (7381.82 u\ + 23673. m? m 2 + 45677.7m? u\ + 31394.4 u x u\ + 

+3819.3 u\) - (169602 m? + 615152 m?m 2 + 1.48093 • 10 6 m?m^ + 

+1.61998 • 10 6 m?m^ + 719990 Mi m^ + 76742 m^), (B12) 

Q 2 = 1- (6mi -2m 2 ) + 

+ (35.5 m? — 7 ' U\u 2 — 30.5 u 2 2 ) + 

-(369.646m? - 130.399 m?m 2 - 306.317mim^ + 164.3 ul) + 
+(7381.82 u\ - 1861.44 m?m 2 - 4739.46 m^ + 2138.26 u x u\ + 
+581.765 u\) - (169602 m? - 33233.6 m?m 2 - 55928.8 m?m^ + 
+53320.3m?m^ - 13607.5mim^ - 6933.24m^). (B13) 



APPENDIX C: RENORMALIZATION-GROUP DIMENSIONS OF THE 
PERTURBATIONS AT THE MN FIXED POINTS 

The analysis of the perturbations at the mn FP is quite simple. In our case m = 2N, 
n = 2 and the relevant symmetry group is 0(2 AT), which is broken by the terms proportional 
to m 2 , w 2 , and w^. If <p a i is the field defined in Eq. (4.6), a — 1, 2, i — 1, . . . 2N, we define 
the following spin-2 and spin-4 operators that transform irreducibly under 0(2iV): 
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K 



(2) 



V aZ,k,l = <t>ai</>aj(/>ak<f>al ~ 



-(pl(5ij(f) ak (j) a i + 5 perm.) 



+ 



1 



4(iV + l)(iV + 2) 
where <p 2 a = £V 0^. Then, the relevant operators are: 

2 1 TV 



2(iV + 2) 

a) 2 (^j4« + S ik 8ji + 5 u 5 jk ), 



Pl — yi ^ yi V aJ+aN,i+pN,j+aN,j+pN' 



a=l a/3=0 ij=l 
1 AT 

^2 = ^ ^ ^l,i+r..V.)+n.\ ! 2 . / 4 I.Y../+ !.V 

a/ 3=0 i?'=l 
1 AT 



K2) T/ 



(2) 



,(2) 



(2) 



a/?7<5=0 ij"=l 



(CI) 
(C2) 

(C3) 
(C4) 
(C5) 



where eoi = — eio = 1 and e o = en = 0. These operators give rise to different breakings of 
0(2iV): 



O(2A0 ■ p(N) <g> 0(2)] © [O(N) ® 0(2)] O(N) <g> [0(2) © 0(2)] 
O(iV) ©S[0(2) ©0(2)]. 
In terms of P 1: P 2 , and P 3 Hamiltonian (1.4) can then be written as 



(C6) 



a a 

+t 3 Pi + UP 2 + t 5 P 3 , 



(C7) 



where 



U\ u 2 1 

(C8) 

Since all operators are irreducible with respect to 0(2N) transformations, if the couplings 
belong to the mn theory, the stability matrix defined with respect to the couplings U has 
the form 

/ Q u Q 21 \ 
tt 21 tt 22 

n= o o fii o o (C9) 

tt 2 
\ tt 2 J 

Note that two eigenvalues are degenerate, since P 2 and P 3 are different components of the 
same irreducible operator V^V^i- The corresponding RG dimensions are Y 1 = — Q 1 , 

Yn = — £lo. 
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